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Summary 


A laminated plate theory analysis is developed to calculate the strain 
energy release rate associated with edge delamination growth in a composite 
laminate. The analysis includes the contribution of residual thermal and 
moisture stresses to the strain energy released. The strain energy release 
rate, G, increased when residual thermal effects were combined with applied 
mechanical strains, but then decreased when increasing moisture content was 
included. A quasi-three-dimensional finite element analysis indicated identi- 
cal trends and demonstrated these same trends for the individual strain energy 
release rate components, Gj and Gjj, associated with interlaminar tension 
and shear. An experimental study indicated that for T300/5208 graphite-epoxy 
composites, the inclusion of residual thermal and moisture stresses did not 
significantly alter the calculation of interlaminar fracture toughness from 
strain energy release rate analysis of edge delamination data taken at room 
temperature, ambient conditions. 
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Nomenclature 


[A] extensional stiffness matrix 

A delaminated area 

a delamination size 

[B] coupling stiffness matrix 

b half width of laminate cross section 

[D] bending stiffness matrix 

E modulus of elasticity 

^LAM longitudinal laminate modulus 

E* longitudinal modulus of laminate completely delaminated along one or 

more interfaces 

Ell »E22 l am ^ na moduli parallel and perpendicular to fiber direction, 
respectively 

F xl ,F x2 ,F x 3 axial f o r c es 

Gj2 lamina shear moduli 

G strain energy release rate 

Gi»Gii»Giii strain energy release rate components due to opening in- 

plane shear, and out-of-plane shear fracture modes 

M m+T m i t i II 

G ,G ,G strain energy release rate due to mechanical, mechanical 

plus thermal, and mechanical plus thermal plus hygroscopic 
loading 

G c critical strain energy release rate for delamination onset 

Gj c critical mode I strain energy release rate for delamination onset 

h ply thickness 

i. length 

N number of plies 

{N} vector of force resultants 

{M} vector of moment resultants 

[Q] transformed reduced stiffness matrix 

t laminate thickness 
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u 

u 

u ek 

u bk 

u ck 

V 

w 

x,y,z 

z k 

“1 

«2 

{ a } k 

M 

*i 

62 


strain energy 
strain energy density 

lamina strain energy density due to extension for the k c ply 

lamina strain energy density due to bending for the k c ^ ply 

lamina strain energy density due to bending-extension coupling for 

the k th ply 

volume 

width 

laminate coordinate system 

distance from laminate midplane to the bottom of the k 1 "* 1 ply 

lamina coefficient of thermal expansion in the fiber direction 

lamina coefficient of thermal expansion normal to the fiber direc- 
tion 

lamina coefficients of thermal expansion in laminate coordinate 
system 

laminate coefficients of thermal expansion 

lamina coefficient of moisture expansion in the fiber direction 

lamina coefficient of moisture expansion normal to the fiber direc- 
tion 


{8> k 


lamina coefficients of moisture expansion in the laminate coordinate 
system 


{6} laminate coefficients of moisture expansion 


AH 

AT 

e 


e ij 

,{e) 


percentage moisture weight gain 

temperature differential from cure temperature to test temperature 
strain 

delamination onset strain 
strain tensor 

vector of in-plane strain components 
mid-plane strain 

strain components due to mechanical, thermal, and hygroscopic 
loading 


3 



vector of curvatures 


M 

{A} 

(u) 

v 12 

°U 

{a} 


laminate coefficients of thermal curvature 
laminate coefficients of moisture curvature 
Poisson's ratio in lamina coordinates 
stress tensor 

vector of in-plane stress components 
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Introduction 


Composite materials are currently being considered for primary structural 
applications that tax their ability to carry load after being damaged. In 
particular, the extensive delamination that may result from low-velocity 
impacts on a thick composite wing skin may severely limit the compressive 
strains the material can sustain, and hence, severely limit the weight saving 
potential of the composite [1,2]. Because of this limitation, there has been 
increased emphasis on developing and measuring improved delamination resis- 
tance in high performance composites [3,4]. One test recently proposed for 
measuring delamination resistance is the edge delaraination tension (EDT) test 
which measures the mixed-mode interlaminar fracture toughness [5,6,7]. This 
EDT test involves measuring the strain at the onset of edge delamination 
during a tensile test of a composite laminate. The measured strain is then 
substituted into an equation for strain energy release rate to obtain the 
interlaminar fracture toughness. Although the mixed-mode nature of the edge 
delamination test has been examined [5, 6, 8, 9], previous analyses have not 
considered the influence of residual thermal and moisture stresses on the 
strain energy release rate as edge delaminations grow. The purpose of this 
investigation was to modify existing analyses to account for the influence of 
residual thermal and moisture stresses on strain energy release rates for edge 
delamination, and to experimentally determine the significance of these ther- 
mal and moisture effects on interlaminar fracture toughness measurement. 
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Strain Energy Release Rate Analysis 
Laminated Plate Theory 

An equation was previously derived for the strain energy release rate, 

G, associated with edge delamination growth in a composite laminate [5]. This 
equation had the following form: 


E^t * 

G - ¥ (E LAM - E > 


0 ) 


where e = nominal tensile strain 

t = laminate thickness 

E L am = longitudinal modulus of the original laminate 

E* = longitudinal modulus of the delaminated laminate 

The G in equation (1) depends on the layup and the location of the delami- 
nated interfaces (which determines E ), but does not depend on delamination 
size. If the undirectional lamina properties Ejj, ^22* v 12» anc * ^12 are 

•fa 

known, then Elam an< * E can ca l cu ^- ate ^ from laminated plate theory and 
the rule of mixtures [5]. The [±30/±30/90/90] g laminate originally proposed 
for measuring interlaminar fracture toughness [5,6,7] delaminates at the 
-30/90 interfaces. After delamination, this layup is modeled as three sublam- 
inates, two [±30 ] 2 laminates and one [ 90 ] 3 laminate, to account for the loss 
in compatibility of transverse contraction as delaminations grow under an 
applied strain. Thus, 

E * = 8E (±30) + 3E (90) (2) 


where E^qj is equal to & 22 ’ an< ^ ^*(±30) can ca l cu l ate d either from 
laminated plate theory or measured from a tensile test of a (±30) ng 
laminate [7]. In the [±35/0/90] g layup, recently proposed as an alternate for 
use in the EDT test [8,9], delaminations grow between 0/90 interfaces. Hence, 
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this layup is modeled as three sublaminates, two [±35/0] g laminates and one 
[90 ] 2 laminate. Thus, 

£ * _ 6E (±35/0) * 2E (90) (3) 

8 

where again E (go) is equal to E 22 » and E (±35/0) may be ca l cu l ate d either 

from laminated plate theory or measured from a tensile test of a (±35/0) s 

laminate. However, assuming the (±35/0) sublaminate is symmetric may lead to 

a small error in E* because of the bending-extension coupling for this 

layup. The axial modulus of an asymmetric layup may be calculated from lami — 

nated plate theory by assuming N , N , k , M , and are all zero 

y xy x y xy 

for a constant e [10,11]. This technique allows for a nonzero ic v and 
x y 

yields a slightly different axial modulus than is calculated for a (±35/0) g 
laminate. This technique was used throughout this study to analyze G. The 
influence of bending-extension coupling on strain energy release rates is 
discussed further in appendix A. 

Previously, delamination onset strains have been substituted into equa- 
tion (1) to calculate the interlaminar fracture toughness. However, equa- 
tion (1) assumes that the delamination growth is due only to the applied 
mechanical strain, e M . Theoretically, however, residual thermal stresses and 
the stresses due to absorbed moisture may contribute to the strain energy 
released as the delamination grows. Therefore, a new analysis was developed 
to account for these thermal and moisture influences. 

Figure 1 shows a simplified model of a composite laminate containing a 
delamination of size "a" growing from either edge. The laminate is assumed to 
be typical of the layups previously used for the edge delamination test, with 
a set of n 90° plies in the center, and some combination of angle plies 
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and/or 0° plies on either side. Delaminations are assumed to occur at the 
interfaces between the central 90° plies and the neighboring angle plies. The 
laminated, LAM, and delaminated regions, SUB and (90) n , are modeled as sepa- 
rate laminates loaded in parallel. The strain energy U is a function of 
delamination size, a, and may be written as the sum of the strain energy of 
three regions: the laminated portion, the separated 90°-ply group, and the 

separated sublaminates on either side of the 90° plies. Hence, 

U(a) = O^^j^(a) + U SUB (a) + Ugg(a) (4) 

For any particular region, the strain energy may be defined in terms of strain 
and stress tensors as 

U = \ I e. .a . .dV (5) 

2 ij ij 

The assumptions of classical laminated plate theory [12] may be used to rede- 
fine the strain energy as 

N h 

u = -T I J- (e}'{cr}dz (6) 

k=l k-1 

where {e}' represents the transpose of the strain vector and N is the 
total number of plies in the laminate or sublaminate. This formulation impli- 
citly accounts for the midplane extension as well as the linear variation in 
stress and strain within each ply arising from out-of-plane bending. For a 
laminate experiencing no bending and having equal ply thicknesses h, 
equation (6) reduces to 

N 

< e >k <‘>k (7) 
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For an asymmetric laminate in uniaxial tension, however, there may be bending- 
extension coupling which requires the inclusion of pure bending and coupling 
terms in the expression for the strain energy. The derivation of these 
bending-extension coupling terms is presented in appendix A. Equation (7) may 
be rewritten as 

U = w£t u (8) 

where u is the strain energy density defined as 


u 


N 


1 

2N 


l 

k=l 


{e}’ {o} k 


(9) 


For the model shown in figure 1, w lam = 2(b-a), w su ^ = 2a, and 

Wg n = 2a. Combining these expressions with equations (4) and (8) yields the 
yu n 

following equations: 


U(a) 2£ [ (t SUB U SUB + C 90 U 90 t LAM U LAM )a + (t LAM U LAM )b J 


( 10 ) 


The strains in equations (6), (7), and (9) are total strains consisting 
of the sura of a mechanical, thermal, and hygroscopic term. Hence, the total 
strain may be written as 

{e} k = {e} k + {e} k + {e} k (11) 
The mechanical component is simply the uniaxial nominal strain, i.e. , 



( 12 ) 
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The residual thermal strains are given by the following: 



( 13 ) 


where a is the coefficient of thermal expansion of the region, 

a is the coefficient of thermal expansion of the k c 1 ply in the 

region, and 

AT is the difference between the stress-free cure temperature and the 
temperature of the ambient environment in which the laminate is 
tested. 

However, for the sublaminates and the (90) n -ply group, ot x becomes a x of 
the original laminate because the delaminated regions are constrained by the 
grips to the original laminates' free thermal expansion.^ 

The hygroscopic (moisture) strains are given by 



(14) 


where 


AH 


is the coefficient of moisture expansion of the region, 
is the coefficient of moisture expansion of the k 1 -* 1 ply in the 
region, and 

is the percent weight difference of the laminate between the time 
of testing and the dry condition. 


^This restriction applies only to the longitudinal strains, e x , and not to 

the transverse, s , and shear strains, 

y 
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However, for the sublaminates and the (90°) n -ply group, (3^ becomes of 

the original laminate because the delaminated regions are constrained by the 
grips to the original laminate free moisture expansion.* 

The laminar stresses in equations (6), (7), and (9) are calculated from 
strains of equation (11) as, 

M k - [qj k {e> k 

_ 1 . 1 . 
where [Q]^ is the transformed reduced stiff qess matrix of the x ply as 

defined in reference 12. 

Recall that for an elastic body containing a planar flaw of area, A, that 
extends under a constant remote strain e, no external work is performed as 
the flaw extends, and hence the strain energy release rate may be written 
simply as 



For the delaminated laminate modeled in figure 1, A = 2£a. Hence, 
dA = 2 l da, and 


JL 

2£ da 


(16) 


Substituting equation (10) into equation (16) and differentiating yields 


G t LAM U LAM t SUB U SUB C 90 U 90 


(17) 


where u for each region is given by equation (9). Therefore, the sprain 
energy release rate analysis requires a ply-by-ply evaluation of the strain 
energy in the laminated and delaminated regions to account for the biaxial 
thermal residual and moisture stresses that may be present in the laminate as 
delaminations grow. 
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Figure 2 shows the influence of residual thermal and moisture stresses on 
the strain energy release rate associated with edge delamination growth in the 
-30°/90° interfaces of the 11 ply [±30/±30/90/90] s laminate. The arrow point- 
ing to the -30/90 interface in the figure caption indicates the delaminated 
interface. Graphite-epoxy lamina properties from reference 13 were used in 
the analysis, along with the following thermal and moisture coefficients [14] 
oj = (-.23 ye/°F) -0.41 yG /°C 
a 2 = (14.9 y e/°F) 26.8 y e/°C 
01 = 0 

8 2 = 5560 pe/weight percent 

A mechanical strain of e = 0.01 and a room temperature condition of 70°F 
were assumed. Therefore, AT was -280°F for the 350°F cured graphite-epoxy 
laminate. The ply thickness was assumed to be 0.0054 inch. As shown on the 
ordinate of figure 2, the strain energy release rate due to mechanical loading 
only, G^, calculated from equation (17) is identical to the value calculated 
from equation (1). If the residual thermal strain is included, the strain 
energy release rate, G^ + ^, is higher than G 1 ^ for the same applied strain. 
However, if the laminate absorbs moisture, the residual thermal stresses are 
relaxed and the strain energy release rate, G^ + ^ + ^, decreases depending on 
the percentage of moisture weight gain, AH. For the case shown in figure 2, 
a moisture weight gain of 0.73 percent completely relaxes the residual thermal 

M4.T4.ll Kf 

stresses, and the strain energy release rate, G , is equal to G where 
only mechanical strains were considered. 

Figure 3 shows a similar relationship for delamination growth in the 
-45/90 interfaces of an 8 ply [+45/0/— 45/90 ] s laminate. As previously shown 
in figure 2, the strain energy release rate for mechanical loading, G**, 
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increases when residual thermal stresses are included, G^ + ^, but decreases 

m+t+h 

when moisture stresses are included, G n 1 . 

Previously [5,6,8,11], the strain energy release rate calculated from 
equation (1) was shown to be equal to the sum of the Gj (interlaminar ten- 
sion) and Gjj (interlaminar shear) components calculated using a quasi- 
three-dimensional finite element analysis [15]. Like the total G, the Gj 
and Gjj components were independent of delamination size; however, the 
percentage of the total G attributed to each component depended on the layup 
[6,8]. In order to determine the accuracy of equation (17) and determine the 
influence of residual thermal and moisture stresses on the individual fracture 
modes, the finite element analysis was modified. 
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Finite Element Analysis 


In the quasi-three-dimensional analysis (Q3D) finite element analysis, 
the laminate is assumed to be long in the x-direction, and every x = constant 
plane deforms exactly in the same fashion as any other x = constant plane. 
Therefore, the gradients of displacements, strains, and stresses, with respect 
to the x coordinate are zero [15-18], and only one x = constant plane must 
be analyzed. Figure 4(a) shows such an x = constant plane. Because of 
symmetries in the problem, only one quadrant, defined by 0 < y < b, 

0 < z < t/2, was analyzed. This region was modeled by eight-noded isoparam- 
etric elements. No singularity elements were used near the delamination 
front. Figure 4(b) shows a finite element model used for this analysis. 

The mesh consisted of 1271 nodes and 386 elements for analysis of the 
[35/— 35/0/90] laminate. The eight-noded elements at the delamination front 
were square, with the dimensions of each side equal to one fourth of a ply 
thickness. Strain energy release rates were calculated using the virtual 
crack extension technique [19]. 

The Q3D delamination problem was solved for three separate loading con- 
ditions: mechanical, thermal, and hygroscopic. For mechanical loading, the 

Q3D analysis is documented in references 16-18. For the thermal and hygro- 
scopic loading the Q3D analysis was performed as presented in appendix B. The 
forces and displacements in the vicinty of the delamination front for the 
combined loading conditions of M+T, and M+T+H, were obtained by appropriate 
superposition of the forces and displacements for the M, T, and H, loading 
conditions. The combined forces and displacements were then used to calculate 
the components of the strain energy release rates Gj, Gjj, and Gjjj. 
sum of three components yielded the total G. 
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Inclusion of the thermal and moisture terms results in a coupled problem 
because of the coupling of the mechanical, thermal, and hygroscopic contribu- 
tions to the strain energy released. Therefore, G results could not be 

2 

normalized by a single loading term such as e , AT, or AH. Hence, the 
results were presented for specific strain levels, at a AT of -280°F, for a 
range of moisture contents, AH. 

Figure 5 compares the finite element results to the laminate plate theory 
results of figure 2 for delaraination in the -30/90 interfaces of the eleven- 
ply [±30/±30/90/90] s laminate. The laminate plate theory results from equa- 
tion (17) agree well with the finite element results. Also shown in figure 5 
is the Gj component of the total G calculated from the finite element 

analysis. A similar trend of the strain energy release rate due to mechanical 
M 

loading, G^ , increasing when residual thermal stresses are included, 

M+t jcf+x+H 

Gj , and then decreasing with increasing moisture content, G^. , is appar- 

ent. The Gj component remained a fixed percentage of the total G for all 
the conditions analyzed. Figure 6 compares the finite element results to the 
laminate plate theory results from figure 3 for delamination in the -45/90 
interfaces of the eight-ply [45/ 0/— 45/90 ] g laminate. Again, the laminate 
plate theory and finite element analysis agree, and the Gj component is a 
fixed percentage of the total G for all the conditions analyzed. However, 
as the following discussion will illustrate, the Gj component may change 
with moisture content for other layups. 

In reference 8, a family of laminates was evaluated where the total 
strain energy release rates for a remote mechanical strain were nearly identi- 
cal, but the Gj percentage varied widely. The three layups examined were 
[45/-45/0/90] g , [45/0/-45/90] s , and [ 0/45/— 45/90] g . A small difference in 
total G was calculated for delaminations in [45/-45/0/90] g and 
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[0/A5/— 45/90] s laminates compared to the [45/0/— 45/90] s laminate because of 

the curvature effect (k^ * 0 ) on sublaminate moduli for the asymmetric 

laminates formed in the first two cases. This curvature effect on moduli is a 

consequence of the non-zero and B 22 terms in the [B] matrix 

influencing the longitudinal modulus. In the third case, the [45/0/-45] 

sublaminate has non-zero B^ and B 25 terms only, corresponding to twist- 

extension coupling. These terms do not affect the longitudinal modulus of the 

sublaminate. Therefore, this sublaminate may be assumed to be symmetric. 

Assuming delaminations occurred between the 90° plies and the adjacent angle 

plies, equation ( 1 ) indicated that was nearly identical for all three 

layups. However, the finite element analysis showed that these three layups 

M 

had relatively high, intermediate, and low percentages, respectively. 

Figure 7 shows the G^ component, normalized by the total G as a function 
of percent moisture weight gain for the three layups. All calculations were 
performed using the delamination onset strains for the appropriate laminates 
as measured in reference [ 8 ]. The solid symbols show the relative ratio of 
Gj to G for the three layups when mechanical strain only is considered. 

If residual thermal stresses are included in the analysis, the Gj/G ratios 
change for the [ 45/— 45/0/90] s and [0/45/-45/90] g layups, but do not change for 
the [45/0/-45/90] layups, as shown by the points labeled M+T on the 
ordinate. Furthermore, the Gj/G ratios continue to change for the first two 
layups as moisture content increases, but G^/G remains unchanged with in- 
creasing AH for the third layup, as shown by the curves labeled M+T+H. 

Figure 8 shows similar results for [ 35/— 35/ 0/90 ] g , [ 35/0/— 35/90 ] s , and 
[ 0/ 35/— 35/90] s laminates assuming a mechanical strain of 0.00254. These 
layups were chosen for the edge delamination test based on a parametric study 
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to minimize the delamination onset strain required to measure the interlaminar 
fracture toughness of a given material {8]. As was the case for the family of 
quasi-isotropic layups studied, the only layup of the three that maintained a 
constant Gj/G ratio when thermal and moisture effects where included was the 
[35/0/— 35/90] s layup that had (0/O/-0) sublaminates after delamination. 
Apparently, layups that form asymmetric sublaminates when they delaminate 
have Gj and Gjj components that vary with AH. For these layups, delami- 
nation is driven by a combination of Poisson mismatch of the individual plies 
and a curvature effect due to bending-extension coupling. However, layups 
that form symmetric sublaminates (or laminates that are asymmetric only in the 
sign of the angle plies) after delamination maintain constant Gj/G ratios as 
moisture content changes. For these laminates, delamination is driven 
entirely by Poisson mismatch. 

Figure 9 shows that the Gj/G ratio also changes with the level of 
mechanical strain. The variation in G^/G with AH is much greater as low 
strain levels than at higher strains. Therefore, the influence of thermal and 
moisture stresses on mixed-mode ratios will not be as significant for tough 
materials that delaminate at relatively high strains as it may be for more 
brittle materials that delaminate at lower strains. 
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Experiments 


Specimen Preparation 

An experimental program was conducted to determine the relative 
significance of including residual thermal and moisture stresses in strain 
energy release rate analysis for interlaminar fracture toughness measure- 
ment. First, a 0.093 m 2 (1 ft 2 ) panel of T300/5208 graphite-epoxy was fabri- 
cated from prepreg tape. The [±35/0/90] layup was chosen to maximize the 
mode I contribution and minimize the possibility of extensive 90°-ply cracking 
occurring before delamination onset [8]. The panel was cured and cooled 
overnight while still under pressure in the vacuum bag in the autoclave. The 
next morning the panel was removed and was immediately cut, using a diamond 
saw, into 44 specimens, 127 mm (5 in.) long by 12.7 mm (0.5 in.) wide.^ Each 
specimen was immediately weighed to the nearest 0.001 gram. Of the 44 speci- 
mens, 4 were tested immediately in the "aftercured condition," 4 were stored 
in a desiccator, 11 were dried in a vacuum oven at 100°C, 8 were soaked in 
water, and the remainder were stored in ambient laboratory air. Specimens 
were chosen from a variety of panel positions for each condition. The per- 
centage moisture weight gain was monitored for each condition for approxi- 
mately 2 months until the specimen weight stabilized. For example, figure 10 
shows the measured moisture weight gain for the specimens that were soaked in 
water. A final specimen weight was measured just before testing. 

Specimens dried in the vacuum oven at 100°C for one month between the 
time they were cut and when they were tested showed an average moisture weight 

^Because the specimens had to be weighed to determine moisture weight gain, 
the specimen dimensions had to be much less than the 254 mm (10 in.) by 
38.1 mm (1.5 in.) size used previously [5-9]. A recent study indicated that 
EDT specimen size did not affect measured delamination onset strains. There- 
fore, an optimum size was chosen that was small enough to weigh but large 
enough to instrument with a 25.4 mm (1 in.) gage length extensoraeter. 
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loss of 0.21 percent. Therefore, an 0.21 percent weight gain was added to 
the AH measured for each condition to approximate the percentage moisture 
weight gain from a truly dry state (table 1). 

Modulus Measurement 

The thickness of each specimen was recorded at three positions along the 
length of each specimen, near either edge, and averaged. Axial strain was 
measured with a 25.4 mm (1.0 in.) extensometer centered over a 76.2 mm 
(3.0 in.) gage length. The blades of the extensometer rested in "V" shaped 
grooves cut into thin aluminum plates that were mounted on the specimen with a 
fast drying adhesive. Load and strain data were recorded on an x-y recorder. 
Table 1 shows the average thickness and modulus measured for the specimens 
tested. Although the average modulus varied somewhat with AH, only the 
water-soaked case showed significant variation. Therefore, a single set of 

*np 

lamina properties were used to calculate and E in equation (1), and 

to calculate strain energy release rates for all but the saturated case. 

The following properties were measured for T300/5208 undirectional lamina 
at room temperature under ambient conditions: 

E u = 18.7 Msi (128 GPa) 

E 22 = 1.23 Msi (8.47 GPa) 

G 12 = 0.832 Msi (5.73 GPa) 
v 12 = 0.292 

An axial laminate modulus of 9.08 Msi (62.5 GPa) was calculated for the 
[ 35/— 35/0/90 ] s layup from laminated plate theory using these properties. This 
calculated modulus was close to the average measured modulus (excluding the 
water-soaked case) of 9.05 Msi (62.4 GPa). The delaminated modulus, E*, 
calculated for delamination in the 0/90 interfaces of the [35/-35/0/90] 

s 

laminate was 7.61 Msi (52.4 GPa). In order to simulate the modulus reduction 
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in the water-soaked case, was reduced to 18.1 Msi (125 GPa). This 

yielded an axial laminate modulus of 8.83 Msi (60.8 GPa) which was the same as 
the measured value. The corresponding value of E for the water-soaked case 
was 7.44 Msi (51.3 GPa). 


Edge Delamination Tests 

For all specimens, an abrupt deviation in the linear load-deflection 
curve ocurred when delaminations formed on the edge. Delamination onset 
strains, e c , for each moisture condition are shown in figure 11. Mean values 
of ply thickness, delamination onset strain, and measured moisture weight gain 
for each condition (shown in table 2) were substituted into equation (17) to 

m+t+h 

calculate G^ . A AT of -280°F was assumed because it represents the 
largest gradient possible for tests conducted at room temperature. 

m+t+h 

Figure 12 shows G^ values plotted as a function of percentage 

moisture weight gain. Also shown in figure 12 are G^ values calculated from 

equation (17) using the average delamination onset strains for each condition. 

M+T+H M 

The G values do not agree with the G values, with differences of 

c c 

the order of thirty percent occurring for the vacuum-dried case (AH = 0%) and 

the water-soaked case (AH = 1.55%). For the ambient condition (AH = 0.52%), 
M+T+H M 

however, the G^ and G^ values were within twelve percent of each other 

because the moisture contribution nearly cancelled the residual thermal con- 


tribution to G. 


In addition to comparing the total G c , the mode I component was plotted 
and compared because previous work had indicated that composites with brittle 
epoxy matrices delaminate under interlaminar tension and are insensitive to 
the presence of interlaminar shear [8,9]. The appropriate mode I percentage 
of the total G c shown in table 2 was determined from the finite element 
analysis for the measured e c and AH (assuming a AT of -280°F), and was 
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used to calculate Gj c (fig. 13). Figure 13 shows the difference between 
Gj c values where thermal and moisture effects have been included and Gj c 

M 

values where they have been neglected. There is less difference between Gj c 
M+T+H 

and G™ values than was observed for the total G„ measurements. At the 

ic c 

ambient conditions (AH = 0.52%), the difference was only seven percent. 

Hence, for most practical purposes, residual thermal and moisture effects may 
be neglected for interlaminar fracture toughness measurements under ambient 
conditions. However, the dry and water-soaked laminates showed significant 
differences in both G c and G lc values using the two analyses (figs. 12 and 
13). Therefore, thermal and moisture effects may need to be included for G 
analysis of composites tested under dry or hot-wet conditions. Furthermore, 
toughened composites containing thermoplastic matrices that are manufactured 
at very high temperatures, but absorb very little moisture, may require 
inclusion of residual thermal and moisture effects to accurately predict G 
values because they may have relatively low moisture contents with relatively 
high AT values. However, as was shown in figure 9, thermal and moisture 
effects become less dominant for toughened-matrix composites that delaminate 
at high strain levels, because the mechanical strain has the greatest con- 
tribution to the strain energy released. 

Discussion 

The significance of residual thermal and moisture stresses to strain 
energy release rates ultimately depends on how these measurements are used. 

If toughness measurements are used to compare materials for improved delamina- 
tion resistance, then these thermal and moisture effects become of secondary 
importance. This is especially true if tests are conducted at room- 
temperature ambient conditions, and the difference in toughness measurements 
for different materials is large [4,5,6]. For example, a seven percent error 
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in G due to neglecting thermal and moisture effects for the [±35/0/90] 

Ic s 

EOT tests is insignificant compared to a ten-fold increase in G^ for com- 
posites with toughened matrices. 

If, however, these interlaminar toughness measurements are used as delami- 
nation failure criteria to predict delamintion growth in composite structures of 
the same material, but with different geometries and loadings, then these 
thermal and moisture effects may become more significant. Other factors may need 
to be addressed to accurately calculate G. For example, the coefficient of 
thermal expansion may not be constant between the cure temperature and room 
temperature. In addition, assuming that the average moisture content of the 
laminate is representative of the moisture content at the delamination front may 
also be in error. Some knowledge of the moisture distribution through the lami- 
nate may be needed. The detailed information required for carefully conducted 
laboratory tests may not be available to analyze the strain energy release rate 
for the delamination growing in the structure. Nevertheless, conducting edge 
delamination tests where these effects can be quantified, and compared to data 
from other interlaminar fracture toughness tests where these effects are not 
present, would help document the relative significance of residual thermal and 
moisture stresses on the interlaminar fracture of composite materials. 

Conclusions 

Based on the analyses and experiments performed during this study, the 
following conclusions have been reached: 

(1) The strain energy release rate calculated for edge delamination 

growth at a constant AT increases when residual thermal stresses 
are included with applied mechanical strains, but decreases when these 
residual stresses are relaxed with increasing moisture content. 
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(2) When residual thermal and moisture stresses are included in strain 
energy release rate analyses, the ratio of Gj to G is constant 
for laminates modeled as symmetric sublaminates after delamination, 
but for laminates modeled as asymmetrical sublaminates after 
delamination the Gj/G ratio changes with moisture content and 
mechanical strain. 

(3) For graphite-epoxy composites residual thermal and moisture stresses 
do not significantly alter the calculation of interlaminar fracture 
toughness from strain energy release rate analysis of edge delamina- 
tion data taken at room temperature and ambient moisture conditions. 
However, for epoxy-matrix composites tested under dry or hot/wet 
conditions the influence of residual thermal and moisture stresses 
may need to be considered to accurately predict strain energy release 
rates for edge delamination. 
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Appendix A 


This appendix presents the formulation of the bending, extension, and 
coupling terms included in the expressions for strain energy of a laminate. 
The strain energy of a given volume may be defined in terms of strain and 
stress tensors as 



£ ij a ij 


dV 


(Al) 


If the volume is composed of orthotropic laminae and is in a state of plane 
stress, the assumptions of laminated plate theory may be used to redefine the 
strain energy in terms of the engineering stress and total strain vectors in 
the laminate coordinate system (x,y,z) as 

N A 

U * ” I / “ {*)' (?) dz (A2) 

k_1 z k-l 

For the general case in which extension and bending occur, the total strain at 
any point is given by the sum of the mid-plane strain and the product of the 
curvature and the distance from the raid-plane 

{e} = { e° } + z{<} (A3) 


The mid-plane strain and curvature are each composed of mechanical, M, ther- 
mal, T, and hygroscopic (or moisture), H, components, and may be written as 


(z°) - (z°) M + { e °) T ♦ |z°) H - (z°} M + («}AT + (s)4H 
M - {«} M + (k) T + (*} H - (k) M + (X)AT + (JJiH (M > 

where {a} and { B } are, respectively, the laminate coefficients of thermal 
and moisture expansion and { A } and {p} are, respectively, the laminate 
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coefficients of thermal and moisture curvature. These vectors are defined as 
follows: 


a) j 'A B- 1 l A B~ _1 (n ) H 

x; AT _B D_ |m) * ly) AH _B D. 1 m) 


-1 /„\ H 


The force and moment resultants, {N} and { M} , are related to the mid-plane 
strain and curvature by 


IK:] C 


where the load and moment are each composed of mechanical, thermal, and hygro- 
scopic components as follows: 


{N} = {N} M + {N} T + {N} H = {N} M + AT l [Q] k {a} k t k 

k=l 


+ AH I I q ] k {e} k t k 


{M} = {M} M + {M} T + {M} H = {m} M + at l [Q] k {a} k (zj - z^J/2 

k=l 

+ 4H l [ ^kl e )t( z k ‘ z k-i ) /2 
k-1 

t* / 

The stress in the k ply is given by the product of the transformed reduced 
stiffness matrix I Q ] k , and mechanical strain of the ply, where the mechanical 
strain is the total strain less the free thermal and moisture strains of that 


Mk = (Ql k (M - {ct} k AT - (B} k AH) 
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here {o}^ and {S}^ are the lamina coefficients of thermal and moisture 
expansion, respectively. 

After substitutions of equations (A3) and (A8), evaluation of 
equation (A2) yields an expression for the strain energy of a laminate in 
terms of the strain energy densities of the individual plies 


U = w£ 


n 

I 

k=l 




ek 


+ u 


ck 


u bk } = 


w l 


N 

I 

k=l 


^“lc 


(A9) 


where 

U ek = 7 { E °MQJ k ({ e °} " M k * T " ^k^ (A10) 

is the strain energy density in the k c ^ ply derived strictly from extension, 

u bk = 6 ( z k + Z kVl + Z k-l) to' (A11) 

is the strain energy density in the k^ ply due to bending only, and 

u ck = i (z k + Vi^M^J 2 ^ 0 } " i°)k AT " Wk AH ) (A12) 

is the strain energy density in the k*-^ 1 ply associated with coupling between 
extension and bending. Hence, for an asymmetric sublaminate, equa- 
tions (A10) - (A12) would be summed in equation (A9), instead of using equa- 
tion (A10) alone, to calculate strain energy density, «SUB» anc * * s use< * 
the strain energy release rate calculation of equation (17). Figure 14 shows 
the influence of including equations (All) and (A12) in the calculation of 

qM+T+H £ or a [-1-35/0/90 ] s laminate delaminated at the 0/90 interface. As 

m+t+h 

moisture content exceeds 0.6 percent, G became significantly larger 

when coupling terms were included. The values tabulated in table 2 were 

calculated with bending and coupling terms included. 
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Appendix B 


This appendix presents the general formulation of the Q3D thermal loading 
analysis. The formulation used for the hygroscopic loading is similar and 
hence is not presented here. 

Consider a laminate subjected to a temperature change of AT and no 
other loading. Due to this temperature change, each ply of the laminate tends 
to either expand or contract, but each ply is constrained by the neighboring 
ply. Therefore, in the interior of the laminate stresses consistant with 
laminated plate theory develop to maintain compatibility of deformation. Near 
free edges, however, interlaminar stresses develop to maintain compatibility. 
Because the laminate is long, every x = constant plane contracts or expands 
in exactly the same manner. However, the amount of this contraction or expan- 
sion is not known. But, because there is no other loading on the laminate, 
the net force in the x-direction must be zero. This condition was used to 
determine the amount of contraction or expansion in the x-direction (and hence 
the strain e^. 

The thermal problem shown in figure 15(a), can be decomposed into two 
subproblems. In the first subproblem (fig. 15(b)) the laminate is subjected 
to a temperature change AT and constrained so that e x = 0. In the other 
subproblem (fig. 15(c)) there is no temperature change, but there is a uniform 
axial strain of e x = 1. 

The subproblem in figure 15(c) is simply the mechanical load problem of 
references 15, 16 and 17. The total solution for the thermal problem in 
figure 15(a) is obtained by scaling the solution for the second subproblem 
figure 15(c) by the unknown axial strain e Q and adding the solutions of the 
two subproblems together. Because the only loading for the problem shown In 
figure 15(a) is the temperature change, AT, the sum of the axial forces of 
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the two subproblems must be equal to zero. This leads to the following 
equation: 


f . + e f 0 = 0 
xl o x2 


which may be solved for the strain due to a temperature change as: 


e o ~ f 


■xl 

x2 


(Bl) 


(B2) 
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TABLE 1.- AVERAGE PHYSICAL PROPERTY MEASUREMENTS FOR T300/5208 

(±35/0/90 ) s LAMINATES 


Environment 

ah% 

From 

Cure 

From 

Dry 

No. of Tests 

t, 

in (mm) 

ElAM’ 

M gi (GPa) 

Vacuum Dried 

-0.210 

0.0 

11 

.0424 

(1.077) 

9.02 

(62.1) 

After Cured 

0.0 

.21 

4 

.0433 

(1.100) 

9.10 

(62.7) 

Dessicator 

.112 

.32 

4 

.0425 

(1.080) 

9.02 

(62.2) 

Ambient 

.310 

.52 

14 

.0427 

(1.085) 

9.08 

(62.6) 

Water Soaked 

1.34 

1.55 

8 

.0428 

(1.087) 

8.83 

(60.9) 
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TABLE 2.- EDGE DELAMINATION DATA REDUCTION FOR CRITICAL STRAIN ENERGY RELEASE RATE 


G c , p S j KJ/M 2 



Environment 

AH% 

h, in 

(mm) 

e c 

g m 

c 

_,Mt-T+H 

U 

C 

G* 

Ic 

„Mt-T+H 

G Ic 


Vacuum Dried 

0.0 

0.0053 (0.135) 

0.00471 

0.674 (0.118) 

1.053 

(0.185) 

0.635 (0.111) 

0.815 (0.143) 


After Cured 

.21 

.0054 ( 

.137) 

.00465 

.670 ( 

.118) 

.962 

( .169) 

.631 ( 

.111) 

.775 ( .136) 


Dessicator 

.32 

.0053 ( 

.135) 

.00493 

.740 ( 

.130) 

.944 

( .166) 

.697 ( 

.122) 

.804 ( .141) 

LO 

U> 

Ambient 

.52 

.0053 ( 

.135) 

.00521 

.825 ( 

.145) 

.937 

( .164) 

.777 ( 

.136) 

.838 ( .147) 


Water Soaked 

1.55 

.0054 ( 

.137) 

.00548 

.866 ( 

.152) 

.628 

( .110) 

.816 ( 

.143) 

.621 ( .109) 







Fig. 3 Influence of residual thermal and moisture stresses 
on energy release rate [45/0/-45/90] T300/5208. 
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Fig. 4 Discretization of Laminate Cross Section for Finite 
Element Analysis 




Fig. 5 Influence of residual thermal and moisture stresses _ 
on strain energy release rate components [+30/+30/90/90] 
T300/5208. t 
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Fig. 6 Influence of residual thermal and moisture stresses on 
strain energy release rate components [45/0/-45/90] 
T300/5208. t b 




Fig. 7 Effect of residual thermal and moisture stresses on percentage 
Mode I strain energy release rate for quasi-isotropic laminate. 
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Fig. 8 Effect of residual thermal and moisture stresses on Mode I 
percentage strain energy release rate. 
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Fig. 9 Effect of Residual Thermal and Moisture Stresses on Gj/G 
as a Function of Moisture Content and Mechanical Strain 













Fig. 14 Strain energy release rate calculation with and 
without bending/extension coupling. 




Fig. 15 Thermal stress analysis for delamination using the 
finite element technique. 
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